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هدف از پژوهش کنونی بررسی دورنمایی از اثرات احتمالی تغییر اقلیم بر جابه‌جایی زمانی تاریخ وقوع 
اولین و آخرین بخبندان‌های پاییزه و بهارهٌ ایران اشعتنت. بدین منظور از داده‌های دیده‌بانی ۳ انششگاه 
همدید کشور (۲۰۱۰-۱۹۸۱) و داده‌های شبیه‌سازی شده ۷۷ علض در دو مدل آب‌وهوای جهانی 
2/21 و ۳13001۷3 در بازه‌های زمانی (۲۰۶۳-۲۰۲۵) و (۲۰۱۸۰۱-۲۰۱۹۹) تحت سه سناریوی انتشار 
۵ ۸2 و 1ظ استفاده کودیت. نتایج. بیانگر جابه‌جایی اولین بخبندان پاییزه به سمت اوایل زمستان و 
جابه‌جایی آخرین یخبندان بهاره به‌سوی اواخر زمستان در گسترة ایران است. پراکنش زمانی-مکانی 
تغییرات متفاوت است؛ به گونه‌ای که بیشترین جابه‌جایی‌های مثبت در رخداد اولین یخبندان پاییزه در دور 
(۲۰۲-۲۰۵) در ایستگاه‌هایی چون حرم‌آباده رشت و گرگان مشاهده می‌شود. میزان تغیبرات در 
ایستگاه‌های شمال شرقی (سبزوار و سمنان)» نیمه جنوبی (کرمان بم و آباده) و بیشتر ایستگاه‌های شمال- 
غرب نسبت به دیگر مناطق کمتتر اسستته: در دوره (۲۰۱۹۹-۲۰۸۰) بیشترین روند منفی در ایستگاه‌های 
گرگان رشت. اردبیل و شهرکرد خواهد بود. خوی» قزوین, بم و کاشان کمترین جابه‌جایی منفی خواهند 
داشت. 


کلیدواژه‌ها: تغییر اقلیم. اولین یخبندان پاییزه» آخرین بخبندان بهاره. مدل ۷۷ 1۵55 ایران. 
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۸۲ جغرافیا و مخاطرات محیطی شماره پیست و چهارم 


۱ مقدمه 

بخش کشاورزی وابسته‌ترین بخش به اقلیم است و اقلیم تعبین‌کنندة اصلی زمان» مکان, منابع تولید و بهره‌وری 
فعالیت‌های کشاورزی است (ریلی, ۹۹۹: ۲۹۸). تاریخ رحداد اولین دمای صفر درجه در پائیز و آحرین رحداد آن 
در بهار به لحاظ کشاورزی دارای اهمیت می‌باشد. این اطلاعات در تعیین گونه‌های مناسب جهت کات در هر 
منطفه‌ای به کار می‌آیند. گرمایش جهانی سبب تغییراتی در شاخص‌های اقلیم کشاورزی مانند: تاریخ رخداد 
یخبندان‌های دیررس بهاره زودرس پاییزه» طول دورهُ یخبندان و طول دورة رشد شده‌است. نتایج بررسی 
پژوهشگران در بسیاری از نقاط جهان در دهه‌های اخیر حاکی از جابه‌جایی زمانی رخداد یخبندان‌ها دارد. کوتر و 
لیدوک! (1۹۹0) با بررسی روند سری‌های زمانی تاریخ‌های رخداد آخرین یخبندان‌های فصل بهار در شمال شرق 
ایالات‌متحدة آمریکا مشاهده نمودند که روند حطی منفی معنادار در تاریخ آغاز یخبندان‌های زودرس فصل بهار دور 
(۱۹۱۱-۱۹۹۰) وجود دارد. مطالعة استرلینگ " (۲۰۰۲) روی روند تاریخ رخداد بخبندان‌ها نشان داد در دورةٌ ۱۹۹4- 
۸ در ایالات‌متحده یخبندان‌های پاییزه دیرتر و یخبندان‌های بهاره زودتر رخ می‌دهد. کیان " و همکاران (۲۰۱۲) به 
این نتیجه رسیدند که زمان وقوع بخبندان‌های دیررس بهاره و زودرس پاییزه در کانادا در طول قرن بیستم به سمت 
فصل زمستان جابجا شده‌اند. ژانگ * و همکاران (۲۰۱6) دریافتند که تاریخ رخداد نخستین یخبندان پاییزه در فلات 
تبت در بازهٌ زمانی (۲۰۱۰-۱۹۲۰) ۱/۵ روز در هر دهه به تعویق افتاده‌است و واپسین یخبندان بهاره ۱/۹ روز در هر 
دهه زودتر رخ می‌دهد. نتیجه بررسی بو "و همکاران (۲۰۱۶) روی تغییرات شاخص‌های بخبندان دورهُ (۱۹۸۰- 
۰ منطقةٌ دریاچه‌های پنجگانة ایالاتمتحده حاکی از روند منفی آحرین بخبندان بهاره و مثبت اولین بخبندان 
پاییزه است که منجر به افزایش طول دورهٌ رشد شده است. پژوهش سکاگز و ایرماک" (۲۰۱۵) حاکی است که فصل 
رشد در ایالات مینسوتا در طول قرن بیستم به سبب جابجایی زمانی در وقوع آخرین یخبندان‌های زودرس و اولین 
یخبندان‌های دیررس افزایش یافته‌است. اما این افزايش در بین ایستگاه‌های مورد مطالعه متفاوت بوده‌است. 

تاریخ رخداد یخبندان‌های پاییزة نواحی جنوبی استرالیا ؛ هفته در دورة (۲۰۱6-۲۰۰۰) نسبت به دهد 1۰ دیرتر 
رخ می‌دهد (کریمپ "و همکاران ۲۰۱7: ۵ امروزه به‌کارگیری خروجی مدل‌های اقلیم جهانی این امکان را به 
وجود آورده تا پژوهشگران دورنمایی از تغییرات احتمالی شاخص‌های آماری یخبندان‌ها را در مقیاس‌های زمانی 


بلندمدت بر بخش‌های مختلف کرةٌ زمین بررسی نمایند. بررسی روند و تغیبرات آيندهٌ شاحص‌های روزهای بخبندان 
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سال ششم شبیه‌سازی ار تغییر اقلیم بر جابه‌جایی زمانی تاریخ وقوع اولین و ... له 


حوضه آبخیز کت اسکیل کوهستان نیویورک ! حاکی از افزايش شاخص‌های دمای کمینه و طول فصل رشد و کاهش 
روزهای یخبندان بود. در دور مشاهداتی (۲۰۰۰-۱۹7۰) در همه حوضه به جز کانون سویل » آخرین بخبندان بهاره 
۳ ۲/۱ روز در هر دهه زودتر و اولین یخبندان پاییزه ۲/۷-۳/۲ روز در هر دهه دیرتر رخ داده و دور رشد 6- 
۶ روز در هر دهه طولانی‌تر شده‌بود. مطابق سناریوهای 0۷1) در دوره‌های ۲۰۶۵-۲۰7۵ و ۲۰۸۰-۲۱۰۰ آخرین 
یخبندان بهاره ۶-۱۱ و ۵-۱۵/ روز در هر دهه زودتر و اولین یخبندان پاییزه ۱-۱۰ و ۶-۱۳ روز در هر دهه دیرتر 
رخ خواهدداد و همچنین دورة رشد ۱۰-۲۵ و ۱۳-۶۰ روز در هر دهمه طولانی‌تر می‌شود (آناندی" و همکاران 
٩4 ۳‏ در حوضهء رودخانه البه" در جمهوری چک در دورءٌ (۲۰۰۰-۱۹۱) و دو دورة آب‌وهوای آینده 
(۲۰۵۰-۲۰۲۱ و ۲۰۷۰-۲۱۰۰) تحت سناریوی ۸13 تاریخ رخداد آخرین یخندان بهاره در دور مشاهداتی ۲ روز 
در هر دهه زودتر به وقوع می‌افتد و اولین یخبندان پاییزه ۱/۸ روز در هر دهه به تأخیر افتاده است. در دهه‌های آتی 
آخرین یخبندان بهاره ب‌طور چشمگیری در مقایسه با آب‌وهوای دورة دیده‌بانی زودتر رخ خواهدداد و اولین یخبندان 
پاییزه به‌طور قابلتوجهی به تأخیر خواهدافتاد (پوتاپ*و همکاران. 7۰۱6 ۳. در ایران نیز بررسی‌های نوحی و 
همکاران (۱۳۸) و طاووسی و درخشی (۱۳۸۸) نشان از جابه‌جایی تاریخ رخداد یخبندان‌های پاییزه و بهاره به 
سمت زمستان دارد. نتایج پژوهش حاصل از شبیه‌سازی‌های مدل‌های گردش عمومی جو برای سال‌های ۲۰۲۵ و 
۰ در ايران نشان داد که تاریخ وقوع اولین یخبندان پاییزه به ترتیب ۵-4 روز و ۸-۱۵ روز به تأخیر خواه دافتاد. 
همچنین رخداد آحرین بخبندان بهاره به میزان 4-۸ روز و ۷-۱۲ روز زودتر رخ می‌دهد و شدت تغییر از شمال به 
جنوب و از غرب به شرق کشور افزایش خواهد یافت (نصیری محلاتی و همکاران ۱۳۸۵: ۲۱. نتایج ارزیابی 
تغییرات یخبندان‌های دیرهنگام خراسان رضوی در دور اقلیمی آینده حاکی از آن است که به‌طور متوسط تاریخ 
وقوع یخبندان‌های دیر رس سبک ۳او یخبندان‌های سنگین ۷ روز زودتر از متوسط تاریخ گذشته به اتمام خواهند 
رسید. همچنین فراوانی وقوع و شدت بخبندان‌های دیررس در سه ایستگاه مورد بررسی کاهش چشمگیری نسبت 
به دوره اقلیمی گذشته خواهندداشت به‌گونه‌ای که میزان احتمال وقوع بخبندان سبک در مشهد ۱۷ درصد 
تربت‌حبدریه ۷ درصد و سبزوار ۳ درصد نسبت به دور گذشته کاهش خواهدیافت این کاهش خطر در مورد 
یخبندان‌های سنگین با 4۷ ۲۳ و ۵۶ درصد به‌ترتیب برای ایستگاه‌های مشهد. تربت‌حیدریه و سبزوار واضح‌تر است 
(اسماعیلی و همکاران. ۱۳۸۹: ۱). شبیه‌سازی رخداد تغییر اقلیم و تأثیر آن بر زمان کاشت. طول دوره رشد و نیاز آبی 


گندم زمستانهٌ بهشهر با استفاده از خروجی مدل‌های ۳3000۷3 و 3 تحت سناریوی 2 در دوره‌های زمانی 
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۲۰۱-۰ ۲۰۱۱-۲۰۵۰ و ۲۰۷۱-۲۱۰۰ حاکی است که با جابه‌جائی تاریخ کاشت از ٩‏ آبان به سمت ۱۹ آذ 
تبخیر و تعرق در تمام دوره‌ها افزایش می‌یابد. علت این افزایش, جابه‌جا شدن مراحل میانی و انتهائی رشد است که 
گیاه در این مرحله حداکثر تبخیر و تعرق را دارد به سمت اواسط بهار که هوا گرم‌تر است بوده و درنتیجه تبخیر و 
تعرق بیشتر خواهدشد (سلیمانی و همکاران. ۱۳۹۱: ۱۱). بررسی تأثیرات تغییر اقلیم بر شماره روز اولین یخبندان 
پاییزه و آحرین یخبندان بهاره در ایران با استفاده از ریز مقیاس نمایی آماری 5125[1 تحت دو سناریوی انتشار ۸2 و 
2 در چند نمونه اقلیمی برای چشم انداز (۲۰۳۹-۲۰۲۰) نشان می‌دهد که میانگین شماره روز آخرین یخبندان بهاره 
در ایستگاه‌های کرمانشاه» گرگان و رشت کاهش و در ایستگاه‌های اصفهان و زاهدان بدون تفاوت و در مابقی 
ایستگاه‌ها افزایش می‌یابد. میانگین شماره روز اولین بخبندان پاییزه در ایستگاه کرمانشاه افزایش, در ایستگاه رشت 
بدون تفاوت و در مابقی ایستگاه‌های مورد مطالعه کاهش خواهد یافت (بهیار و همکاران ۱۳۹۲: ۱۱۷). پژوهش‌هایی 
که در ایران صورت گرفته با یک یا چند ایستگاه در منطقه‌ای محدود صورت گرفتهاست و اگر هم سراسر ایران را 
پوشش داده با یک مدل و سناریو و تعداد معدودی ایستگاه انجام یافته‌است. در این پژوهش سعی می‌شود انرات 
افزایش دمای کمینه بر جابه‌جایی زمانی تاریخ وقوع اولین و آخرین بخبندان‌های پاییزه و بهارة ایران در دهه‌های آتی 


با دو مدل و سه سناریو شبیه‌سازی شود. 


۲ منطقه مورد مطالعه 
کیلومتر مربع است. میانگین سالانه دمای ایران حدود ۸ درجه سانتی گراد است. اگرچه در دهه‌های اخیر میانگین دما 


اتایشن باه این یراتفر کیفوه آنران به یک شک و اندازه ره ورین فان افدانی اکن در آیران 
در دمای کمینه سك داده‌است (علیجانی و همکاران. 2۱۳۹۰ ۳-۲ 


۳ مواد و روش‌ها 

داده‌های مورد نیاز پژوهش شامل دو گروه داده‌های مشاهداتی و داده‌های شبیه‌سازی شده‌است. داده‌های 
مشاهداتی شامل: مقادیر بارش, تابش, دمای بیشینه و کمينة روزانةٌ ۶۳ ایستگاه همدید در نواحی گوناگون ایران است 
که دارای داد کامل سی سال (۲۰۱۰-۱۹۸۱) بودند» از سازمان هواشناسی کشور دریافت شده است (شکل ۱). گروه 
دوم داده‌های شبیه‌سازی شلده دورة آینده (دهة۲۰۵۰ و ۲۰۸۰) با استفاده از برونداد ریزگردانی شده دو مدل 
آب‌وهوای جهانی 0۳07۷21 و 11200۷13 در نسخه ۵ مدل 1۸126-1۷0 تحت سناریوهای انتشار ( ,۸2 ,131 
8 تولید شدند. 
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شکل ۱. پراکنش ایستگاه‌های مورد مطالعه ایران 


مدل 1۸135-۷۷0 برای مدل‌سازی متغیرهای هواشناسی توزیع‌های آماری پیچیده‌ای به کار می‌برد. سری‌های 
فوریه درجه حرارت را تخمین می‌زنند. دمای کمینه و پیشینهٌ روزانه به صورت فرایندهای تصادفی و با میانگین و 
انحراف معیارهای روزانه که وابسته به وضعیت تر یا خشک بودن روز مورد نظر هستند. مدل‌سازی می‌شوند. سری 
فورية مرتبهٌ سوم برای شبیه‌سازی میانگین و انحراف معیار دمای فصلی به کار می‌رود. رفتار حاکم بر مقادیر مانده‌ها که 
از تفریق مقادیر میانگین از مقدار دیده‌بانی شده به دست می‌آیند. از طریق تحلیل خود همبستگی زمانی داده‌های کمینه 
و بیشینه به دست می‌آیند (سمنوف و بارو ۲۰۰۲ 6). مدل در سه مرح واسنجی» ارزیابی و ایجاد داده‌های 
هواشناسی انجام می‌شود. نیاز اساسی در مرحلهٌ واسنجی, فایلی است که مشسخص‌کننده رفتار اقلیم در دورهُ گذشته 
می‌باشد. این فایل با استفاده از داده‌های مشاهداتی دورة پایه تهیه شده و مدل بر اساس آن اجرا می‌شود. برای 
مدل‌سازی آینده در ابتدا مدل بایستی آگاهی و شنات کافی از رفتار آماری دوره دیدبانی را داشته باشد. این شنانعت 
با واردکردن داده‌های روزانةٌ دیده‌بانی به مدل در بخش ارزیابی صورت می‌گیرد که نهایتاً منجر به تولید فایلی با 
پسوند ۷8 می‌شود. در بخش ارزیابی داده‌های آب‌وهوایی دیده‌بانی شده برای تعیین ویژگی‌های آماری‌شان تحلیل 
می‌شوند. درصورتی که نتایج به دست آمده از مقايسة داده‌های تولیدی مدل با داده‌های دیده‌بانی قابل قبول باشد؛ می- 


توان آن‌ها را در سناریوهای آب‌وهوایی آینده مورد استفاده قرار داد. برای ارزیابی این‌که توزیع احتمال داده‌های تولید- 
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۸ جغرافیا و مخاطرات محیطی شماره پیست و چهارم 


شده نزدیک به توزیع احتمال داده‌های دیده‌بانی ایستگاه‌های مورد بررسی است يا نه؛ توزیع احتمالاتی داده‌های دیده- 
بانی و داده‌های شبیه‌سازی‌شده در گزين 150 0 با استفاده از آزمون نکویی برازش کلم وگرف اسمیرنف و میانگین 
و انحراف معیار داده‌ها با آزمون تی نمونه‌ای زوجی (تی وابسته) بررسی و معناداری آن‌ها در سطح معناداری ۰/۰۵ 

در این پژوهش برای بررسی عدم قطعیت از روش بوت‌استراپ استفاده شد. روش بوت‌استراپ برای برآورد 
میزان دقت برآوردگرها؛ هنگام مشاهدات مستقل کاربرد دارد (ایران پناه و محمدزاده ۱۳۸۶: 1۵۳). برای محاسبة 
فاصله اطمینان به روش بوت‌استراپ در نرم‌افزار 5۳55 ابتدا سری زمانی ماهانة دمای کمينة ۳۰ ساله ایستگاه‌ها 
(۲۰۱۰-۱۹۸۱) با استفاده از حروجی (02501106) لارس تهیه شد و محدودهة بالا و پایین معناداری میانگین و انحراف 
معیار دمای کمينةٌ ماهانٌ هر ایستگاه در سطح ۸۹۹ محاسبه و با ميانگین و انحراف معیار دمای کمن داده‌های دیده‌بانی 
مقایسه شدند. هرگاه مقادیر برآوردی مدل‌ها در محدود؛ اطمینان داده‌های مشاهده‌ای قرار گیرد بیانگر اطمینان در 
سطح اعتماد مورد نظر می‌باشد و درصورتی که خارج از محدوده باشد. بیانگر عدم اطمینان به مقدار برآوردی است. 

به دلیل وجود عدم قطعیت در پیش‌بینی‌های اقلیمی ناشی از تفاوت‌های ساختاری در مدل‌های اقلیمی جهانی و 
همچنین عدم قطعیت در تغییرات شرایط اولیه یا پارامتری کردن مدل استفاده از چند مدل و سناریو در پیش‌بینی‌ها 
توصیه می‌شود (سمنوف و ستارتونوویج؛ ۲۰۱۰ ۲) اجرای چند مدل و سناریوی گسترةٌ وسیعی را برای تحلیل ایجاد 
می‌کند. با انجام این کار می‌توان عدم قطعیت را در تولید داده‌های اقلیمی آینده به حداقل رساند. 

در پیش‌فرض مدل ۱۵1۸1۹5-۷۷6 مدل با سناریوهای انتشار حانوادهءٌ ۸و ظ وجود دارد. لازم است جهعت 
انتخاب کاراترین مدل‌هاء ابتدا عملکرد و کارایی این مدل‌ها را بر مبنای سناریوهای انتشار برای شبیه‌سازی داده‌های 
آینده مورد ارزیابی و تأیید قرارگیرد. انتخاب براساس اجرای مدل در دوره‌های متفاوت با بالاترین ضریب تعیین (*8) 
و پایین‌ترین شاخحص خطا شامل: حطای جذر میانگین مربعات (18۷15۳) و شاحص آماری توافق ویلموت (4) که 


مقداری بین صفر تا یک است که مقدار یک بیانگر بهترین برازش است. صورت گرفت (روابط ۱ تا ۳). 
بت . , 1 
(مم(نهم یج( 
(۲۱-۲) 2 (ته)ر << 


در این رابطه مقدار 2 پیانگر ارتباط خطی بین داده‌های شبیه‌سازی‌شده و مشاهداتی بوده که مقدار ان بین صفر ۳ 


1 
طٌ 


تانطق ( 


یک می‌باشد. هرچه این مقدار به یک نزدیک‌تر باشد نشان‌دهنده رابطه حطی قوی‌تر بین دو مقدار می‌باشد. 20و لا 


نکر این اه سای و شیتازی که تورسظ ال کاو ک نانگم کل دیا اس 
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1۳ 
۱ ۲ 1 
رابطة۲ (تعطامتصنع) 2 ]۳۷5 
ای 9 
ددم تیصنو) ,2 
رابطةٌ ۳ سا 0-1 
(تعطم-ت عطم| + وطم تصنو|) ,ر 


در این روابط 51170,1 مقادیر شبیه‌سازی شده 008 مقادیر مشاهداتی و 01981 میانگین مقادیر مشاهداتی است. 
عدم کارایی مدل در برآورد مقادیر را نشان می‌دهد. در شاحص آماری توافق ویلموت () هم مقادیر بین صفر تا یک 
است که مقدار یک بیانگر بهترین برازش است (هاشمی. ۱۳۹۶: .)٩۲‏ از مجموعٌ ۱۵ مدل لارس, دو مدل آب‌وهوای 
دما دارا بودند. انتخاب شدند. 

۳ ۱ مدل۳۲۱۵۸۱۷۲3 

مدل 11۸190۷13 از سال ۳۲۰ روزه استفاده می‌کند و دارای یک شبکه فضایی با ابعاد ۲/۵ درجه عرض 
جغرافیایی در ۳/۷۵ درجه طول جغرافیایی می‌باشد. این مدل به عنوان یکی از کامل‌ترین و مشهورترین مدل‌های 
گردش عمومی جو در نظر گرفته می‌شود. مدل منحصربه‌فرد است و در آن نیازی به تنظیمات شار جهت تولید یک 
سناریوی واقعی نیست (گوردون" و همکاران ۲۰۰۰ ۱۶۷). 

۳ ۲. مدل ۶۳0۷۲21) 

این مدل از مجموعهٌ مدل‌های آزمایشگاه دینامیک شارهای ژئوفیزیک دانشگاه پرینستون آمریکاست. مدل برای 
شبیه‌سازی آب‌وهوا جوی و اقیانوسی در مقیاس زمانی روزانه طراحی شده‌است. مدل اتمسفری نسخهٌ 003۷21 
و عرض جغرافیایی. دارای ۵۰ تراز عمودی است (دلورث " و همکاران ۲۰۰7: 14۶). 

داده‌های شبیه‌سازی شده با این دو مدل برای بازه‌های زمانی (۲۰۳۵-۲۰۶) و (۲۰۹۹-۲۰۸۰) تحت سه 
سناریوی انتشار ۵2 ۸113 و 131 اجرا و تولید شدند. بعد از تشکیل بانک اطلاعاتی داده‌های مشاهداتی و شبیه‌سازی 
شده بازة زمانی ۳۰۳ روزه جهت مطالعه انتخاب شد؛ یعنی اول شهریور به عنوان روز مبنا انتخاب و روزهای دیگر به 


ترتیب نسبت به این مبداً شمارش و آخرین روز ۳۱ خرداد در نظر گرفته شد. پس از تنظیم داده‌ها بر اساس روز 


1 0 
2 ۷0 


۸ جغرافیا و مخاطرات محیطی شماره پیست و چهارم 


ژولیوسی تاریخ رخداد نخستین یخبندان پاییزه و وایسین یخبندان بهاره برای هر ایستگاه به صورت جداگانه 
استخراج گردید. سپس برای مقايسة تغییرات احتمالی تاریخ رخداد اولین و آخرین یخبندان در دهه‌های آتی اخحتلاف 
داده‌های مشاهداتی و شبیه‌سازی در اقلیم میانی و دور با دو مدل و تحت سه سناریو محاسبه و نقشههای تغییرات 


تاریخ رخداد اولین و آخرین بخبندان در بازه‌های زمانی مورد نظر ايران بر مبنای داده‌های ایستگاهی تهیه ک 0 


6 بحث و نتایج 

بررسی ننایج روایی مدل لارس در شبیه‌سازی دادده‌های کمينة دما با آزمون 5 بیانگر این است که احتلاف معنی 
داری بین داده‌های شبیه‌سازی شده و دیده‌بانی وجود ندارد و با آزمون تیء اختلافات به غیر از ایستگاه بجنورد در ماه 
ژانویه ایستگاه‌های بوشهر و شهرکرد در ماه مارس. ایستگاه اصفهان در ماه فوریه و ایستگاه زاهدان در ماه اوت در 
سطح ۰/۰۵ معنی‌دار نیست. 

در کل می‌توان گفت که مدل)181*5-۷۷ از توانایی قابل قبولی جهت شبیه‌سازی داده‌های دمای کمينة ایستگاه- 
های مورد مطالعه در ايران برخوردار است و خطای موجود تصادفی است. 

تحلیل عدم قطعیت میانگین دمای کمینه ایستگاه نمونهُ سنندج در شکل (۲) آورده شده‌است. نتایح بیانگر آن است 
که ميانگین دمای کمينة ایستگاه سنندج ۱۰ ماه و انحراف معیار دمای کمينه ۲ ماه در محدودة اطمینان ۸۹٩‏ قرار دارد. 
در کل ایستگاه‌ها میانگین دمای کمینه در ۷۲/۷ درصد ماه‌هاه ولی ۷ درصد انحراف معیار ماه‌ها در محدوده اطمینان 
قرار دارد. درحالی‌که در تمام ایستگاه‌ها بیشتر ماه‌ها انحراف معیار داده‌های دیده‌بانی فاصلة زیادی با سطح اطمینان 
داده‌های برآوردی ندارد. نتایج حاکی از ضعف مدل در برآورد انحراف معیار است که برای کاهش عدم قطعیت 


آن‌چنان‌که ذکر شد از دو مدل و سه سناریو بهره گرفته شده‌است. 


5و0 - - - 0.99 600۲ ۲ 


میانگین کمینه دما - ساسیوس 
۲ 
6 مه 


۱۵۷ 00۲ 5۴۴ عنام الا لالاا ۱۸۸۷ ۸۴۶ ۱۸۴ ۴۶5 لا۸۵( 


شکل 9 نمودار بوت استراپ میانگین و انحراف معیار دمای کمينة ایستگاه سنندج در محدودة اطمینان ۹۹ 
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بررسی نتایج روند افزايش میانگین سالائة تغییرات دمای کمینه و بيشينة ایران در دهه‌های ۲۰ ۵۰ و ٩۰‏ آتی با 
میانگین مدل‌های ۲1۵00۷13 و 00721 نسبت دور (۲۰۱۰-۱۹۸۱) نشان می‌دهد که میزان افزایش دمای کمینه 
بیشتر از افزایش دمای بيشینه خواهد بود جدول (۱). افزايش دمای کمینه. تأثیر زیادی روی شاخص‌های یخبندان 
چون: فراوانی روزهای یخبندان؛ تغییر در تاریخ رخداد اولین و آخرین یخبندان‌ها» طول دورة یخبندان و طول دور 


رشد دارد. 


جدول ۱. میانگین سالانةٌ افزایش دمای کمینه و بيشينة ایستگاه‌های مورد بررسی ایران در دهه‌های آتی با میانگین 


مدل‌های 1120613 و 610121 و سناریوهای انتشار ۸18 ,۸2 و191 نسبت دور (۲۰۱۰-۱۹۸۱) 


دوره سناریو دمای بیشینه 7 دمای کمینه ۶) 

۸۷۱ ۷ ۸1 

۷۵ ۳۱ ۸72 ۲۰۱۱-۳۰ 
۰/۳۷6۹ ۰/۵ ۳1 

۱۳/۳۷ ۱۳/۳ ۸11 

۲۳/3 ۳/۳ 2 ۲۰۵ 
۱/۸ ۱/۳۷5۹ 121 

۳/۵ ۳/۷ ۹۱۰ 

1۹ 1۷ ۳ ۲۰۸۰-۲۵ 
۱۳/۳۷ ۱۳/۳ 1۳1 


میانگین رخداد اولین روز یخبندان در ایران 

نقشه پهنه‌بندی میانگین اولین روز یخبندان در ايران نشان می‌دهد که میانگین اولین روز بخبندان در حد فاصل ۲۰ 
مهر تا ؛ آبان در ارتفاعات شمال غربی (سقز. همدان اردبیل و زنجان» شمال شرقی (بجنورد تربت‌حبدریه و 
بیرجند و همچنین در ارتفاعات زاگرس مرکزی (شهرکرد) رخ می‌دهند و دیرترین میانگین آن نیز در حد فاصل ۲۰ 
دی تا ؛ بهمن کمی دورتر از سواحل جنوبی ایران در نوار باریکی موازی با ساحل (از ایرانشهر تا دزفول) و قسمت- 
هایی از سواحل شمالی (از بابلسر تا بندر انزلی) رخ می‌دهد. بیشتر بخش‌های سواحل جنوبی ایران بدون یخبندان 
هستند؛ بنابراین فاصله بین زودترین و دیرترین میانگین اولین روز بخبندان در ایران ۱۰7 روز می‌باشد (شکل ۳. نتایج 


۹۰ جغرافیا و مخاطرات محیطی شماره بیست و چهارم 


اولین روز یخبندان (۲۰۱۰ 1۹۸۱ 


۵آبان - ۱۹ آبان 


۰ مهر - ۴ آبان 


44 46 48 50 52 54 56 58 60 62 84 


شکل ۳. میانگین رخداد اولین روز یخبندان در ايران در بازهٌ زمانی (۲۰۱۰-۱۹۸۱) 

در اقلیم میانی (۲۰3۵-۲۰۶) بر اساس مدل 00۳6121 و سناریوهای انتشار ۸13 ,۸2 و131 در ایستگاه‌هایی 
مانند گرگان» خرم‌آباده رشت و اردبیل بیشترین تغییر نسبت به دورة پايهةٌ (۲۰۱۰-۱۹۸۱) دیده می‌شود. سراسر نوار 
ساحلی جنوب با فاصلهٌ زیادی از ساحل از چا‌بهار و ایرانشهر در جنوب شرق تا آبادان و دزفول در جنوب غرب 
بدون یخبندان خواهد بود. ميانگین تغییرات به ترتیب سناریوهای یادشده به ترتیب ۸۰۱۳و ۷ روز خواهد بود. در 
مدل 112061۷13 میانگین تغییرات ٩‏ ۸ و ۷ روز خواهد بود. روند مثبت در تمام ایستگاه‌ها به غیر از سقز دیده می- 
شود. در سقز در سناریوی 131 دو روز نسبت به دوره پایه زودتر رخ خواهد داد. درحالی‌که براساس مطالعة بهیار و 
همکاران (۱۳۹۲) اولین یخبندان پاییزه پیش‌بینی‌شده در ایستگاه‌های کرمان زاهدان. شیراز و گرگان زودتر رخ 
خواهدداد. 

در دورهُ (۲۰۹۹-۲۰۸۰) بر اساس مدل 0700۷121 و سناریوهای ۸13 ,۸2 و۳81 ۲۱۰۱۸ و ٩‏ روز و در مدل 
3 (۱ ۱ و ٩‏ روز روند مثبت نسبت به دور (۲۰۱۰-۱۹۸۱) دیده می‌شود. تغییرات مکانی تا حدود زیادی 
همانند دورة قبلی است فقط فراوانی آن بیشتر شده است. بیشترین جابه‌جایی‌ها در ایستگاه‌هایی چون خرم‌آباد رشت 
و گرگان مشاهده می‌شود. میزان تغییرات در ایستگاه‌های شمال شرقی سبزوار و سمنان و نیمه جنوبی کرمان بم و 
آباده و بیشتر ایستگاه‌های شمال‌غرب نسبت به دیگر مناطق کمتر است شکل (4). در پژوهش بهیار و همکاران 


(۱۳۹۲) نیز ایستگاه تبریز تغییر زیادی نداشته‌است. 


سال ششم شبیه‌سازی اثر تغییر اقلیم بر جابه‌جایی زمانی تاریخ وقوع اولین و ... ٩‏ 


۲۰۱۵ 
۳1 


0۳۵221 


۱۳2۱9 ۵/6 


۲۷ ۰ ۴ 


۲۰۵-۸۰ 


"0۳/21 


۱۰۳-96 


5۲ 1 


۴ ۱ ۱ لا تا وا سا اه ۲ 


یدون 
یخبندان 


شکل 1 تغییرات نخستین روز بخبندان در ایران در دوره‌های (۲۰۳۹۱۵-۲۰۶۲) و (۲۰۹۹-۲۰۸۰) نسبت به دوره 
مشاهداتی (۲۰۱۰-۱۹۸۱) براساس خروجی مدل‌های 6۳21 و 112460113 تحت سه سناریوی انتشار ۸11 


2 و 1ظ1. 


۹۲ جغرافیا و مخاطرات محیطی شماره بیست و چهارم 


میانگین رخداد آخرین روز بخبندان در ایران 
زودترین ميانگین آحرین روز بخبندان در ایران در نوار باریکی کمی دورتر از سواحل شمالی خلیج فارس و 
دریای عمان و موازی با سواحل در اوایل بهمن رخ می‌دهد. همچنین دیرترین روز میانگین سللانهٌ پایان یخبندان در 


محمودی و همکاران (۱۳۹۲) در این زمینه تا حدود زیادی منطبق با نتایج پژوهش فعلی است. 
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شکل ۵. ميانگین رخداد آخرین روز بخبندان در ایران در بازة زمانی (۲۰۱۰-۱۹۸۱) 
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در بازةٌ زمانی (۲۰۵-۲۰۶۲) روند کاهشی یعنی جابه‌جایی به سمت زمستان در همه ایستگاه‌ها دیده می‌شود. بر 
اساس مدل ۳0۷21) و سناریوهای انتشار ۸1 ,۸2 و81 آخرین بخبندان بهاره با میانگین ۱7 ۱۵ و ٩‏ روز و در 
مدل 1120013 ۱۶,۱۳ و ۱۱ روز زودتر رخ خواهدداد. در ایستگاه‌های گرگان رامسر بابلسر و رشت بیشترین 
جابه‌جایی رخ خواهدداد. در ایستگاه‌های بم و زابل و تبریز میزان جابه‌جایی کم خواهدبود. نکتة قابل توجه این‌که 
ایستگاه‌های غربی چون همدان سقز سنندج و زنجان در رخداد اولین بخبندان پاییزه دارای تغییر کمی بودند در 
وقوع آخرین بخبندان بهاره دارای جابه‌جایی زیادی هستند؛ یعنی میزان جابه‌جایی رخداد آخرین پخبندان بهاره به 
سمت زمستان در آن‌ها بیشتر از جابه‌جایی اولین یخبندان پاییزه به سمت زمستان است. 

در دهةّ ۲۰۹۹-۲۰۸۰(۸۰) میزان تغییر بر اساس سناریوهای مذکور در مدل 0۳0۷21 به ترتیب ۲۰ ۲۲ و ۱۳ و 
در مدل ۲1200۷۲3 ۲۰,۲۷ و ۱1 روز خواهد بود. بیشترین روند منفی در ایستگاه‌های گرگان» رشت. اردییل و 
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شکل 1 تخییرات واپسین روز بخبندان در ایران در دوره‌های (۲۰۳۵-۲۰۶) و (۲۰۱۹۹-۲۰۸۰) نسبت به دوره 
مشاهداتی (۲۰۱۰۱-۱۹۸۱) پراساس خروجی مدل‌های 01 و ۲۱20۷3 تحت سه سناریوی انتشار ۸۵۸18 


2 و 1ظ1. 


۹ جغرافیا و مخاطرات محیطی شماره پیست و چهارم 


۵ نتیجه گیری 

نتایج مشاهده شده در این مطالعه با بسیاری از یافته‌های پیشین در خصوص تغییرات تاربخ آغاز و پایان یخبندان‌های بهاره 
و پاییزه در ايران و سایر نقاط جهان هم‌خوانی دارد و بیانگر جابه‌جایی رخداد یخبندان‌های پاییزه و بهاره به سمت زمستان 
است. به گونه‌ای که میانگین رخداد اولین یخبندان پاییزه ایران در دورة (۲۰۵-۲۰۶) بر اساس مدل 1:01۷]21) سناریوهای 
انتشار ۵۸113 ۸2 و 131 به ترتیب ۸۰۱۳و ۷ روز و در مدل ۲1۵0۵03 ٩‏ ۸و ۷روز دیرتر رخ خواهد داد. در دورهٌ (۲۰۸۰- 
۵۹ بر اساس مدل ۳0۷21 و سناریوهای ۸113 ,۸2 و1 ۱۸ ۲۱ و ٩‏ روز و در مدل ۳1200713 ۱۰ ۲۱ و ٩‏ روز 
روند مثبت نسبت به دور (۲۰۱۰-۱۹۸۱) دیده می‌شود. وایسین یخبندان بهاره در اقلیم میانی بر اساس مدل 2۳۷121) و 
سناریوهای انتشار ۸18 ,۸2 و91 آخرین بخبندان بهاره با میانگین ۰۱1 ۱۵ و ٩‏ روز و در مدل 112007713 ۱۶,۱۳ و ۱۱ روز 
زودتر رخ خواهد داد. در دههٌ ٩۰‏ (۲۰۹۹-۲۰۸۰) میزان تغییر بر اساس سناریوهای مذکور در مدل ۳0۷]21) به ترتیب ۲۰ 
۲و ۱۳ و در مدل 1200713 ۲۰۰۲۷ و ۱ روز خواهد بود. یافته‌های پژوهش تفاوت‌هایی از نظر جابه‌جایی زمانی و مکانی 
با پژوهش‌های پیشین دارد. در مطالعة نصیری محلاتی و همکاران (۱۳۸۵) در اقلیم میانی (۲۰۵۰) تاریخ وقوع اولین بخبندان 
پاییزه در ایران ۸-۱۵ روز به تأخیر خواهدافتاد. همچنین رخداد آخرین یخبندان بهاره ۷-۱۲ روز زودتر رخ می‌دهد. 
پراساس یافته‌های این پژوهش در اقلیم میانی روند مثبت در تما ایستگاه‌ها به غیر از سقز دیده می‌شود. درحالی‌که 
ترآساشن مطالم هار و هکاران (۱۳۹۲) وین بان بابوه بترشتی ند دز ایستگاه‌های کرسان اهاز نارای 
گرگان زودتر رخ خواهدداد. تفاوت‌هایی که وجود دارد ناشی از انتخاب نوع مدل. سناریو احتلاف در انتخاب طول 
دورهُ مشاهداتی و شبیه‌سازی است. 

این جابه‌جایی زمانی از سویی منجر به کاهش دورة یخبندان و کاهش تولید گندم دیم می‌شود. از سوی دیگر با افزایش 
طول دورهةٌ رشد باعث افزایش تبخیر و تعرق» تخلیةٌ رطوبت خاک و کاهش جریان منجر به تغییر چرخه هیدرولوژیک می‌شود 
و به‌علاوه باعث نیاز آبی بیشتر گیاهان و درنهایت استفادة بیشتر از آب می‌شود. این تغییرات و اثرات زیان‌بار آن می‌تواند 
تهدیدی برای دستیابی به توسعه پایدار در کشورمان باشد. ازاین‌رو افزایش آگاهی از چگونگی تغییرات در پراکندگی زمانی_- 
مکانن: شاعصی‌های نان هر افعاد اس اشی‌های مار کارا اه عم وهی فر کنو و اسرآن کته پر دو‌مرهن 
پیامدهای مخرب آن است نیازمند مطالعه و توجه بسیار است. 
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